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Abstract.  We  propose  a  method  for  treating  the  Dirichlet  boundary  condi¬ 
tions  in  the  framework  of  the  Generalized  Finite  Element  Method  (GFEM). 
We  are  especially  interested  in  boundary  data  with  low  regularity  (possibly 
a  distribution).  We  use  approximate  Dirichlet  boundary  conditions  as  in  [11] 
and  polynomial  approximations  of  the  boundary.  Our  sequence  of  GFEM- 
spaces  considered,  Id  =  1,2,...  is  such  that  Sm  (X  and  hence  it 

does  not  conform  to  one  of  the  basic  FEM  conditions.  Let  h M  be  the  typi¬ 
cal  size  of  the  elements  defining  S ^  and  let  u  E  Hrn~^1(D)  be  the  solution  of 
the  Poisson  problem  —  Aw  =  /  in  Q,  u  =  0  on  dD,  on  a  smooth,  bounded 
domain  Q.  Assume  that  \\v^\\u1/2{dQ.)  —  ff1(Q)  f°r  ^ 

and  | u  —  wi|#i(Q)  <  Ch™\\u\\Hm+i(n),  u  E  Hrn+1(D)  fl  Hq  (Q),  for  a  suit¬ 
able  ui  £  Sp.  Then  we  prove  that  we  obtain  quasi-optimal  rates  of  con¬ 
vergence  for  the  sequence  E  S M  of  GFEM  approximations  of  u ,  that  is, 
\\U  ~  Un\\ <  Ch™\\f\\Hrn-l(Qy  We  also  extend  our  results  to  the  inho¬ 
mogeneous  Dirichlet  boundary  value  problem  —A u  =  f  in  f2,  u  =  g  on  <9f2, 
including  the  case  when  f  —  0  and  g  has  low  regularity  ( i.e is  a  distribu¬ 
tion).  Finally,  we  indicate  an  effective  technique  for  constructing  sequences 
of  GFEM-spaces  satisfying  our  conditions  using  polynomial  approximations  of 
the  boundary. 
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Introduction 

In  the  past  few  years,  meshless  methods  for  the  approximation  of  solutions  of 
partial  differential  equations  have  received  increasing  attention,  especially  in  the 
Engineering  and  Physics  communities.  The  reasons  behind  the  development  of 
such  methods  are  the  difficulties  associated  to  the  mesh  generation,  particularly 
when  the  geometry  of  the  domain  is  complicated.  As  in  the  case  of  the  usual  Finite 
Element  Method,  one  of  the  major  problems  in  the  implementation  of  meshless 
methods  is  the  enforcement  of  Dirichlet  boundary  conditions.  It  is  the  purpose  of 
this  paper  to  address  the  problem  of  enforcing  Dirichlet  boundary  conditions  in  the 
Generalized  Finite  Element  Method  framework.  We  are  especially  interested  in  the 
case  when  the  Dirichlet  boundary  data  has  low  regularity,  including  the  case  when 
it  is  a  distribution  on  the  boundary. 

The  classical  Rayleigh-Ritz  methods  for  elliptic  Dirichlet  boundary  value  prob¬ 
lems  assume  that  the  trial  subspace  functions  fulfill  the  boundary  conditions.  Nev¬ 
ertheless,  the  construction  of  such  subspaces  implies  many  difficulties  in  practice 
when  the  boundary  of  the  domain  is  curved.  Therefore,  several  approaches  are 
known  for  dealing  with  the  Dirichlet  boundary  conditions  on  domains  with  curved 
boundaries.  One  approach  is  to  modify  the  variational  principles  by  adding  ap¬ 
propriate  boundary  terms  so  that  there  will  be  no  need  for  the  trial  subspaces  to 
fulfill  any  condition  at  the  boundary.  See  the  works  of  Babuska  [2,  4],  Bramble  and 
Nitsche  [12],  and  Bramble  and  Schatz  [13,  14],  among  others,  for  examples  of  how 
this  approach  works  in  practice.  Another  approach  (used  also  in  this  paper)  is  to 
use  subspaces  with  nearly  zero  boundary  conditions.  This  ideea  was  first  outlined 
by  Nitsche  [27]  and  further  studied  by  Berger,  Scott,  and  Strang  [11]  and  Nitsche 
[28]. 

Yet  another  approach  is  the  Isoparametric  Finite  Element  Method  or  IFEM  with 
curved  finite  elements  along  the  boundary.  See  [18]  and  references  therein,  or  [17], 
[19,  21,  23,  24,  32,  33],  among  many  others,  for  more  recent  work  and  applications. 
This  approach  is  typically  used  in  connection  with  a  numerical  quadrature  scheme 
computing  the  coefficients  of  the  resulting  linear  systems.  In  the  applications  of  this 
method,  except  in  special  cases  (such  as  when  O  is  a  polyhedral  domain)  the  interior 
O h  of  the  union  of  the  finite  elements  is  not  equal  to  fi,  although  the  boundary  of 
0 h  is  very  close  to  Oil.  That  is,  the  approximate  solution  Uh  is  sought  in  a  subspace 
Vh  C  Hq( ilh)  and  so,  the  homogeneous  Dirichlet  boundary  condition  u  =  0  on  Oil 
is  “approximated”  by  the  boundary  condition  Uh  =  0  on  dilh-  In  fact,  iq,  is  the 
solution  of  a  variational  equation  a,h(uh,Vh)  =  ( fh,Vh)h  for  all  Vh  €  Vh,  where 
a/l(-,-)  is  a  bilinear  form  which  approximates  the  usual  bilinear  form  defined  over 
H1(ilh)xH1{i^h),  and  fh  €  Vf*  approximates  the  linear  form  Vh  €  14  — >  Jn  fvhdx, 
where  /  is  an  extension  of  /  to  the  set  il h  ■ 

Our  approach  has  certain  points  in  common  with  the  isoparametric  method  just 
mentioned  in  the  fact  that  we  are  using  polynomial  approximations  of  the  boundary. 
However,  our  method  does  not  require  non-linear  changes  of  coordinates.  Our 
method  combines  the  approaches  in  the  papers  of  Berger,  Scott,  and  Strang  [11] 
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and  Nitsche  [28].  Our  definition  of  the  discrete  solution  is  as  in  [11],  whereas  our 
assumptions  are  closer  to  those  of  [28] .  We  have  tried  to  keep  our  assumptions  at  a 
minimum.  This  is  possible  using  partitions  of  unity,  more  precisely  the  Generalized 
Finite  Element  Method  or  GFEM,  a  method  that  originated  in  the  work  of  Babuska, 
Caloz,  and  Osborn  [6]  and  further  developed  in  [4,  7,  8,  10,  22,  25,  26,  34], 

Our  construction  is  different  from  the  IFEM  in  that  we  do  not  require  compli¬ 
cated  non-linear  changes  of  coordinates.  Moreover,  our  method  uses  non-conforming 
subspaces  of  functions  and  it  does  not  have  to  deal  with  extensions  over  larger  do¬ 
mains.  It  is  closely  related  to  [9]  which  uses  GFEM  for  elliptic  Neumann  boundary 
value  problems  with  distributional  boundary  data.  The  GFEM  is  a  generalization 
of  the  meshless  methods  which  use  the  idea  of  partition  of  unity.  This  method 
allows  a  great  flexibility  in  constructing  the  trial  spaces,  permits  inclusion  of  a 
priori  knowledge  about  the  differential  equation  in  the  trial  spaces,  and  gives  the 
option  of  constructing  trial  spaces  of  any  desired  regularity.  We  mention  that  the 
GFEM  is  also  known  and  used  under  other  names,  such  as:  the  method  of  “clouds,” 
the  method  of  “finite  spheres,”  the  “X-finite  element  method,”  and  others.  See  the 
survey  by  Babuska,  Banerjee,  and  Osborn  [4]  for  further  information  and  references. 

Let  us  now  describe  the  main  results  of  this  paper  in  some  detail.  Let  ft  C  R" 
be  a  smooth,  bounded  domain  with  boundary  dft.  Let  /  €  H~1(ft)  and  u  £  H1(fi) 
be  the  unique  solution  of  the  Poisson  problem 

(1)  —Ait  =  f  on  ft,  u  =  0  on  dft. 

(For  most  of  our  results,  we  shall  assume  that  /  €  Hp~1(fl)  with  p  >  1.) 

Assume  that  we  are  given  a  sequence  hlt  — >  0  and  a  sequence  C  H1(fl)  of  test- 
trial  spaces.  The  parameters  h M  play  the  role  of  the  size  of  the  elements  defining 
S^.  We  define  the  discrete  solution  £  S M  in  the  usual  way:  Blu^,  v^)  =  (/,  uM) 
for  all  v M  £  Sp,  where  B(u,v )  :=  (Vu,  (see  Equation  (3)  below).  We  do 

not  assume  S M  to  satisfy  exactly  the  Dirichlet  boundary  conditions,  that  is,  we  do 
not  assume  C  H^ft). 

Let  us  fix  from  now  on  a  natural  number  m  £  N  =  {1,2,...}  that  will  play,  in 
what  follows,  the  role  of  the  expected  order  of  approximation.  We  shall  make  the 
following  two  basic  assumptions.  The  first  assumption  is  that  our  approximating 
functions  satisfy  Dirichlet  boundary  conditions  approximately: 

•  Assumption  1,  nearly  zero  boundary  values:  ||iv|| H1/2(dn)  fs  Cl/i™||iyl||jji(n) 
for  any  £  S^. 

The  second  assumption  is  an  approximation  assumption  that  will  be  used  also  for 
non-homogeneous  boundary  conditions.  For  that  purpose,  let  us  consider  a  second 
sequence  of  subspaces  S M  C  H1(ft),  S M  C  <SM,  which  are  not  required  to  satisfy  any 
boundary  conditions. 

•  Assumption  2,  approximability:  for  any  u  £  HJ+1(ft),  0  <  j  <  m,  then 
there  exists  it/  £  S f  such  that  |it  —  «■/ li/qn)  <  Ch^\\u\\jjj+i^Qy  If  u  =  0  on 
d ft,  then  we  can  take  ui  £  S 

These  two  assumptions  are  formulated  in  more  detail  in  Section  1. 

Our  paper  is  divided  into  two  parts.  In  the  first  part,  consisting  of  the  first 
four  sections,  we  prove  some  general  approximation  results  for  the  Poisson  problem 
with  Dirichlet  boundary  conditions.  The  approximations  (or  discrete  solutions)  be¬ 
long  to  some  abstract  spaces  (for  zero  boundary  conditions)  or  S M  (for  general 
boundary  conditions)  that  are  required  to  satisfy  certain  reasonable  assumptions 
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(Assumptions  1  and  2,  for  zero  Dirichlet  boundary  conditions,  and  Assumptions  1, 
2,  and  3  for  non-zero  boundary  conditions).  In  the  second  part  of  the  paper  we 
provide  examples  of  Generalized  Finite  Element  Spaces  satisfying  these  assump¬ 
tions.  Under  Assumptions  1  and  2,  our  main  approximation  result  in  Section  1  is 
the  following. 

Theorem  0.1.  Let  S M  C  H1( O)  be  a  sequence  of  finite  dimensional  subspaces 
satisfying  Assumptions  1  and  2  for  a  sequence  — »  0  and  0  <  p  <  m.  Then 

the  (unique)  solutions  u  and  u M  of  Equations  (1)  and  (3),  respectively,  with  f  £ 
ffp-1(£2)  satisfy 

\\u  ~  M/JffRfi)  -  Hp+1(Q)  <  ChP\\f\\Hp-i^, 

for  constants  independent  of  p  and  f  £  Hp~1( 0). 

In  Section  2  we  extend  our  results  to  the  non-homogeneous  Dirichlet  boundary 
conditions  case  u  =  g  on  OQ,  with  g  £  Hm+1//2(dfl).  In  essence,  we  will  be  looking 
for  a  sequence  Gk  of  approximate  extensions  of  g ,  that  is,  a  sequence  of  elements  of 
Hm+1(Ll)  subject  to  the  following  assumption.  Recall  that  the  sequence  /iM  should 
be  thought  of  as  the  “typical  size”  of  the  elements  defining  Sj,  and  satisfies  hfl  —>  0. 
•  Assumption  3,  approximate  extensions:  There  exists  a  constant  C  >  0 
such  that,  for  any  g  £  Hm+1/2(dfl),  there  exists  a  sequence  Gk  £  Sk 
such  that  ||Gfc|an  —  g\\Hl/2(dn)  —  Ch™\\g\\Hm+i/2(Qri)  an(l  ll^fcll-H^+qn)  < 

C||fl,ll/r"*+V2(an)- 

Let  Wk  be  the  exact  solution  of  —Awk  =  f  +  A  Gk  in  O,  Wk  =  0  on  <90.  Also, 
let  ( Wk)/j.  G  Sji  be  the  discrete  solution  of  this  equation,  namely,  the  solution  of  the 
discrete  variational  problem 

(2)  B((wk)n,v)  =  (/  +  AGfc,u)i2(fl),  v  £  S^, 

where  /  £  Hm~1( ft)  is  the  data  of  Equation  (13).  The  result  we  prove  in  Section  2 
is  the  following. 

Theorem  0.2.  Suppose  Assumptions  1,  2,  and  3  are  satisfied.  Let  Uk  '■=  ( wk)k  + 
Gk-  Then  there  exists  a  constant  C  >  0  such  that  the  solution  u  £  Hm+1(Ll)  of 
Equation  (13)  satisfies 

llw  -  uk\\m(n)  <  C'/i™(||/||ffm-i(o)  +  ||5llffm+i/2(an)). 

In  order  to  deal  with  low  regularity  boundary  data,  in  Section  3  we  consider 
the  Dirichlet  problem  —  Aw  =  /  in  fl,  u  =  g  on  dfl,  with  g  £  H1/2~s(dfl )  and 
/  £  s  >  0.  Thus  both  /  and  g  may  be  distributions.  We  say  that  u  = 

(uq,u\)  £  Hl~s(Ll)  :=  1L1_S( fl)  ®  H~1/2~s{dLl)  is  a  weak  solution  of  the  problem 
—Att  =  /  in  fl,  u  =  g  on  the  boundary  <90  if  B(u,v)  =  —{f,v)n  +  (g,duv)dci,  for 
all  v  £  H1+S(Q).  Then,  the  result  of  Section  3  is  the  following. 

Theorem  0.3.  Let  g  £  7L1/2_S(90)  and  f  £  7L_1_S(0).  Then  there  exists  a  unique 
weak  solution  u  =  (uo,Ui)  £  iJ1_s(0)  for  the  problem  —A u  =  f  in  O,  u  =  g  on 
the  boundary  dLl  and  this  solution  satisfies 

||uo||ffi-»(n)  +  llMill ff-l/2-qan)  <  Cn,s(||fl||_H-i/2-s(aa)  +  ||/||fl--i-»(n))i 
for  a  constant  Cq  s  that  depends  only  on  O  and  s. 
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The  second  part  of  this  paper  is  dedicated  to  constructing  concrete  examples  of 
Generalized  Finite  Element  Spaces  S M  and  satisfying  the  Assumptions  1,  2,  and 
3  of  the  first  part.  In  fact,  we  will  prove  that  these  assumptions  are  easy  to  fulfill 
with  a  “flat-top”  partition  of  unity  and  polynomial  local  approximation  spaces.  The 
exact  conditions  are  formulated  in  Section  4.  The  proof  that  the  resulting  GFEM 
spaces  satisfy  the  Assumptions  1,  2,  and  3  is  in  Section  5.  In  addition,  in  Section  6 
we  also  prove  interior  estimates  for  the  error  u  —  Uk,  where  u  is  the  solution  of 
the  distributional  boundary  value  problem  —A u  =  0  in  fi,  u  =  g  on  dfl,  with 
distributional  data  g  £  H1^2~s(dfl),  s  >  0,  and  u £  S &  is  the  discrete  solution  in 
our  GFEM  spaces.  The  last  section  contains  some  comments  and  a  discussion  of 
some  further  problems. 

In  this  paper,  we  shall  use  the  convention  that  C  >  0  indicates  a  generic  constant, 
independent  of  g,  which  may  be  different  each  time  when  used,  but  is  independent 
of  the  free  variables  of  the  formulas. 

Part  1.  Approximate  Dirichlet  boundary  conditions 

1.  Homogeneous  boundary  conditions 

In  this  section,  we  give  a  proof  of  Theorem  0.1.  We  begin  by  fixing  the  notation 
and  then  we  prove  some  preliminary  results.  In  particular,  /  £  H~1( £2)  and  u  £ 
Hg(fl)  is  the  solution  of  the  Poisson  problem  (1). 

Recall  that  f l  C  is  a  smooth,  bounded  domain ,  fixed  throughout  this  paper. 
We  shall  fix  in  what  follows  m  £  N  =  {1,2, . . .},  which  will  play  the  role  of  the 
order  of  approximation.  We  want  to  approximate  u  with  functions  £  S^,  g  £  N, 
where  S M  C  H1( Q)  is  a  sequence  of  finite  dimensional  subspaces  that  satisfy  the 
Assumption  1  and  2  formulated  next.  In  those  assumptions,  the  sequence  h^  — >  0 
should  be  thought  of  as  the  “typical  size”  of  the  elements  defining  £M.  Our  first 
assumption  is: 

•  Assumption  1  (nearly  zero  boundary  values).  There  exists  C  >  0 
such  that 

11^11^1/2(00)  <  C'/i”1||uM lltfi(fj)  for  any  v M  £  S„. 

So  S M  does  not  necessarily  consist  of  functions  satisfying  the  Dirichlet  bound¬ 
ary  conditions.  Let  |u|ffi(o)  :=  [Jn  \S/u\2dx]1^2 .  Our  second  assumption  will  also 
be  used  for  non-homogeneous  boundary  conditions,  so  we  formulate  it  also  for  a 
sequence  of  spaces  S M  C  Hl{VL),  S M  C  Stl . 

•  Assumption  2  (approximability):  There  exists  C  >  0  such  that  for  any 
0  <  j  <  to,  any  u  £  FP+1(fl),  and  any  g  £  N,  there  exists  uj  £  such 
that 

If  u  =  0  on  dil,  then  we  can  take  ui  £  S 

We  now  proceed  to  the  proof  of  Theorem  0.1.  We  first  need  some  lemmas.  We 
begin  with  the  following  classical  result  [1,  16]. 

Lemma  1.1.  For  v  £  H1( fi)  there  is  a  constant  C  that  depends  only  on  S2  such 
that 

IMIffqa)  ^  C[MfU(n)  +  llw  lli,2(3n)]  • 
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From  this  lemma  we  obtain  that  and  ||wM||/j-i(n)  are  equivalent  norms 

on  with  equivalence  bounds  independent  of  p. 

Lemma  1.2.  There  exists  C  >  0  such  that 

C  1kiu|Ff1(n)  <  ll^vllffqn)  ^  ^\vr\ h1^) 
for  all  p  large  enough  and  all  G  S^. 

Proof.  From  Lemma  1.1,  we  have 

IWIffbn)  ^  C[\vr\H1(n)  +  II*vIIl2(c>q)]  S  + 

where  the  last  inequality  is  a  consequence  of  Assumption  1.  Therefore,  for  p  large, 
h n  is  small  enough  and  we  get 

ll^||_f/i(n)  <  C{  1  -  Ch2™)~1/2 |uM|ffi(n) 

which  is  enough  to  complete  the  proof.  □ 


Lemma  1.2  allows  us  now  to  introduce  the  discrete  solution  u „  of  Equation  (1) 
using  the  standard  procedure.  Let  B(v,w)  :=  Jn  Vv  ■  Vwdx  be  the  usual  bilinear 
form.  For  p  large,  let  us  define  the  discrete  solution  G  S M  of  the  Poisson  problem 
(1)  by  the  usual  formula 

(3)  =  /  f(x)vl_l(x)dx,  for  all  v p  G 

Jn 

Let  v  be  the  outer  unit  normal  to  dfl  and  dS  denote  the  surface  measure  on  dfl. 
Similarly,  let  G  S^,  for  p  large,  be  the  solution  of  the  variational  problem 


(4)  B(wll,vli)  =  /  dI/u(x)vlj,(x)dS(x),  for  all  G  S^, 

Jon 

where  u  is  the  solution  of  Equation  (1).  Note  that  we  need  Lemma  1.2  to  justify 
the  existence  and  uniqueness  of  the  (weak)  solutions  u M  and 

Lemma  1.3.  Let  u  be  the  solution  of  the  Poisson  problem  (1)  and  let  and  w M 
be  as  in  Equations  (3)  and  (4).  Then  B(u  —  u^  —  w^^Vy)  =  0  for  all  G  hence 

I u  -  Uf,  -  <  \u  -  for  all  v M  G  S„. 

Proof.  This  is  obtained  from  Assumption  1  as  follows 


(5)  B(u,Vp)  =  /  Vu-Vv^dx  =  /  fv^dx  +  /  (d„u)vildS(x)  =  B(u^ 
J  £1  J f2  J  dd 


+  «V,vM), 


for  all  Vfj,  G  S^. 


□ 


We  now  proceed  to  estimate  u M  and  wfl . 

Lemma  1.4.  Let  u  be  the  solution  of  the  Poisson  problem  (1)  and  let  w M  be  the 
solution  of  the  weak  problem  (4).  Then,  for  p  large,  we  have 

(6)  1 1  II  h1  (n)  <  Ch™\\u\\Hi(Q), 

with  C  a  constant  independent  of  p  and  u. 
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Proof.  Let  us  first  assume  that  /  £  L2(12).  Then  u  £  H2(Q)  and  hence  d„u  is 
defined  and  belongs  to  P[1^2{dkl).  Green’s  formula  gives  jgnduu{x)w{x)dS{x)  = 
/n  (Attic  +  Vm  •  \7w)dx.  Since  the  map  H1( 12)  — >  H1^2(d 12)  is  surjective  with  a 
continuous  splitting,  we  see  that 

(7)  <  C,|lwll/r1(n)- 

By  continuity,  this  inequality  will  continue  to  hold  for  f  £  iJ_1( 12). 

We  have 

\\wij-\\‘h1(Q)  <  C\wp  1^1(0)  =  CB(wp,Wp)  =  C  dvu{x)wp{x)dS{x) 

J  dQ, 

<  C\\dvu\\H-i/2(dnj)\\'w^\\Hi/2(d£i)  <  G/i™||Tt||ffi(n)||Tc/i||ffi(n)- 
Therefore  ||tcM||_H-i(n)  <  Ch™\\u\\Hi(n),  as  claimed.  □ 

From  this  lemma  we  obtain  the  following. 

Lemma  1.5.  For  /i  large,  the  solution  u p  of  the  weak  problem  (3)  satisfies 

(8)  IWIffhn)  <  C,||w||H1(n)i 

with  C  a  constant  independent  of  p  and  u. 

Proof.  Let  us  first  observe  that  Lemmas  1.2  and  1.3  and  Equation  (7)  give 

ll^vll/rqn)  A  C'\uh\h1(Q)  =  CB(up,  Up)  =  C\B(u,  Up)  —  Bfwp,  tt^)] 

=  C[B(u,u„)  -  (<9„u,  Up)dn\  <  C[\B(u,Up)\  +  \{duu,u^)on\] 

<  G||tt||t/i(o)||um||tti(o)  +  C\\dvu\\H-i/^^air)\\up\\Hx/2^dQ) 

<  C,|MU1(n)IWI.H'1(n)  +  C,^™IIuIIh1(^)IIumIIh1(^)- 

Now  it  is  easy  to  see  that  11^11^1(0)  <  G||Tt||/ri(o),  as  claimed.  □ 

We  are  ready  now  to  prove  Theorem  0.1. 

Proof.  We  shall  assume  p  =  m ,  for  simplicity.  The  proof  in  general  is  exactly  the 
same.  Lemma  1.3  and  the  projection  property,  together  with  Lemma  1.4,  give 


(9)  | u  —  u^\h^(q)  <  \u  —  Up  —  Wp\Hi(n)  +  |«vl h1^) 

<  | u  —  u/|/fi(Q)  +  Chnp\\u\\Hi(Q)  <  G/i(7||u||ffm+i(o), 

where  for  the  last  line  we  also  used  the  approximation  property  (Assumption  2). 

The  estimate  in  the  H1- norm  is  obtained  from  Lemma  1.1,  Equation  (9),  As¬ 
sumption  1,  and  Lemma  1.4  as  follows 

||w  -  Up  llifi(fi)  <  C\\u  -  Up  Iffi(fi)  +  1 1  Up  1 1 7,2  (50)  ] 

<  Ch^\\u\\Hm+i^  +  Ch™ ||?vl|ffi(n)  <  Ch™\\u\\Hm+i(ny 

The  proof  is  now  complete.  □ 


In  view  of  some  further  applications,  we  now  include  an  error  estimate  in  a 
“negative  order”  Sobolev  norm.  We  let  H~l( ft)  to  be  the  dual  of  Hl( ft)  with  pivot 
L2(12).  Since  12  is  a  smooth  domain,  can  also  be  described  as  the  closure 

of  C°°  (12)  in  the  norm 


IMIff-i(fi)  =  sup 
0^0 


|(n,  4>)l2( 0)1 


(10) 
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(Note  that,  in  several  other  papers,  H  l(fl)  denotes  the  dual  of  Hl0{ 17).) 

Theorem  1.6.  Let  0  <  l  <  in,  0  <  p  <  m,  and  7  =  min{7  +  p  +  1 ,  m}.  Then, 
under  the  assumptions  of  Theorem  0.1,  the  solutions  u  and  u M  of  Equation  (1)  and 
Equation  (3),  respectively,  satisfy 

llM  -  <  Ch^WuWHP+HQ)  <  ^11/11^-1(0), 

for  a  constant  C  >  0  independent  of  p  and  f  €  Hp~1(  17). 


Proof.  The  proof  of  this  theorem  is  an  adaptation  of  the  usual  Nitsche-Aubin 
trick.  Indeed,  let  us  denote  by  F  £  ff1+2(fl)  the  unique  solution  of  the  equation 
—A F  =  (f,  F  =  0  on  <917,  for  (f>  £  arbitrary,  non-zero.  Then  there  exists 

a  constant  C  >  O,  independent  of  <j>,  such  that  ||F|| ff(+2(Q)  <  C\\<j)\\Hi(ny  By 
Assumption  2,  there  exists  Fj  £  S/  such  that 

(11)  \F  -  Fi\m(n)  <  Ch1^1  ||F|| Hl+2(Q)- 

Then,  the  inequality  (11)  leads  to  the  following  easy  observation,  which  will  be  used 
later, 

(12)  |A/|ffi(o)  =  \F—(F— .F/)|Hi(Q)  <  |A|ffi(o)  +  C'/i|I+1||F||/i-i+2(o)  <  C\\4>\\hi^q). 

In  the  following  calculation,  we  shall  use  Equation  (5)  in  the  first  inequality,  and 
then  Theorem  0.1,  Equations  (11)  and  (12),  1.4,  and  1.5  for  the  last  inequality,  to 
obtain 


II  |(u  -  UM,</>)L2(n)| 

W  -  Un\\H-‘(n)  =  sup - TT-rr; - =  sup 

<^o  Il0l|j?'(fi) 


<  sup 


\B(u-  u^,F  -  F^l 


B(u  -  Up,  F)  +  fdQ  u^dvFdS 
Ion  u/j.9vFdS  | 


</>Jo  ll'/llz/qo)  <t>4 a  ll^llsHn) 


<^o  Il0l|ff‘(n) 

\u  —  Ufi,\Hirm\F  —  E/|ffi(Q)  |w/i|i?i(n)|A|//i(n) 

<  sup - Arr - -  +  sup - iAA - — 

0  lAllfl-qa)  <t>^  0  \\<P\ \w(n) 


sup 

<t>*  0 


KmIU2(3»)II  AE||L2(3n) 

IAIA!(T2) 


<  C'/i^+!+1||u||j?p+i(n)  +  C'/i”"||u|Ai(n)  +  C,/i™llulli/i(n) 


<  Ch^WuWHp+nQ), 


by  the  definition  of  7.  This  completes  the  proof. 


□ 


For  p  =  0,  the  proof  of  the  above  result  requires  the  full  strength  of  Assumption 
1.  The  case  p  =  0  is  the  one  needed  for  the  results  of  Section  6. 


2.  Non-homogeneous  boundary  conditions 


In  this  subsection  we  provide  an  approach  to  the  non-homogeneous  Dirichlet 
boundary  conditions.  That  is,  we  consider  the  boundary  value  problem 


(13) 


—Aw  =  /  on  17, 

u  =  g  on  <917. 


Our  approach  is  to  reduce  it  to  the  case  g  =  0  and  then  to  use  the  results  on  the 
Poisson  problem  (1).  In  a  purely  theoretical  framework,  this  is  achieved  using  an 
extension  G  of  g  and  then  solving  the  problem  —Am;  =  /  +  AG,  w  =  0  on  <917. 
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The  solution  of  (13)  will  then  be  u  =  w  +  G.  This  gives  that  the  problem  (13)  has 
a  unique  solution  u  £  Hp+1( ft)  for  any  /  £  Hp~1( ft)  and  g  £  H1^2+P(d ft)  and  it 
satisfies 

IMI/fp+qn)  <  C'(||/||ffp-i(Q)  +  ||p||ffi/2+P(c)n)), 
for  a  constant  C  >  0  that  depends  only  on  f t  and  p  £  Z+.  (This  result  is  valid  also 
for  p  =  0.) 

In  practice,  however,  we  need  to  slightly  modify  this  approach  since  it  is  not 
practical  to  construct  the  extension  G  (this  is  especially  a  problem  if  g  has  low 
regularity,  that  is,  if  g  is  a  distribution,  for  instance).  We  will  be  looking  therefore 
for  a  sequence  Gk  of  approximate  extensions  of  g,  that  is,  satisfying  the  following 
assumption.  Recall  that  the  sequence  hfl  should  be  thought  of  as  the  “typical  size” 
of  the  elements  defining  and  satisfies  h;l  — >  0. 

Assumption  3  (approximate  extensions).  We  assume  that  there  exist  a 
sequence  of  spaces  Sk,  Sk  C  Sk,  satisfying  Assumption  2  and  a  constant  C  >  0 
such  that,  for  any  g  £  Hm+1/2(dfl),  there  exists  a  sequence  Gk  £  Sk  such  that 
||Gfc|an  —  #||/fi/2(3fi)  <  Ch'fl\\g\\Hm+i/2(dn')  and  ||Gfc||//™+i(n)  <  G||<7||ffm+i/2(an). 

The  proof  of  Theorem  0.2  follows  below. 

Proof.  Remember  that  Wk  was  introduced  as  the  exact  solution  to  the  boundary 
value  problem  —A  Wk  =  f  +  A  Gk  in  ft,  Wk  =  0  on  dft.  Let  (wk)f_i  £  S M  be  the 
approximate  solution  of  this  equation,  as  in  Equation  (2). 

We  have  that  Vk  :=  Wk  +  Gk  solves  the  boundary  value  problem 

—Avk  =  /  on  f l,  Vk  =  Gk  on  df t. 

Hence  the  difference  u  —  Vk  solve  the  boundary  value  problem  A (u  —  Vk)  =  0, 
( u  —  Vk)  =  g  —  Gk  on  df t.  From  this  and  Assumption  3  we  obtain 

(14)  \\u  —  VftHfl-i(n)  <  C\\g  -  Gk\\H1/2(dQ)  ^  Ch1k\\g\\Hm+i/2(gny 
Theorem  0.1  and  Assumption  3  then  give 

ll^fc  —  =  || Wk  —  (wfc)fe||/fi(Q)  <  Ch™\\f  +  AGfc||//m-i(n) 

<  Ch™  (\\f\\Hm-l(p,)  +  ||Gfc||Hm  +  l(Q))  <  G/l”(||/||ffm-l(n)  +  llplljym+J/a^n)). 
Hence 

(15)  ||ufe  —  uk ||ffi(n)  =  \\wk  -  («Jfc)fc||i?i(n)  <  G/i^(||/||Hm-i(Q)  +  ||flilflm+i/2(an))- 

From  Equations  (14)  and  (15)  we  obtain  ||it  —  WfcUnqn)  <  G/i™(||/||j/m-i(Q)  + 
2 (an)))  which  is  what  we  had  to  prove.  □ 

3.  Distributional  boundary  data  and  the  “inf-sup”  condition 

Let  us  consider  the  Dirichlet  problem  (13)  ( i.e —A u  =  f  in  ft  and  u  =  g  on 
dfl,)  with  g  £  Hl/2~s (dft)  and  /  £  H~1~s(fl),  s  £  R.  If  s  <  0,  it  is  well  known  that 
the  boundary  value  problem  (13)  has  a  unique  solution  u  £  H1~s(fl).  Moreover, 
there  is  a  constant  GqjS,  depending  only  on  ft  and  s  <  0,  such  that 

I|w||h1-s(q)  <  GniS(||/||#-i-S(Q)  +  ||3||_f/i/2-»(ao)). 

In  this  section,  we  extend  the  above  result  to  the  case  when  g  £  H1//2~s(dfl), 
with  s  >  0.  Our  approach  is  based  on  the  so  called  “inf-sup”  condition  [3]. 
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3.1.  The  weak  formulation.  Let  us  define  the  functional  space  := 

fL1_s( Cl)  ®  H~1/2~a{dCl).  Intuitively,  for  an  element  u  =  (u0,ui)  G  H1~s( fi), 
the  first  component  u$  should  be  thought  of  as  u  in  the  interior  of  Cl,  while 
the  second  component  u±  represents  the  normal  derivative  dvu  on  Oil.  Let  B  : 
Hl~s(Cl)  x  H1+s(Cl)  — >  C  be  the  bilinear  functional  defined  by 

B(u,v)  :=  (u0,Av)q  +  («i,v)an, 
where  u  =  (u0,Ui)  G  £f1_s( fl). 

Definition  3.1.  Let  g  G  H1^2~s(dCl)  and  /  G  H~1~s(Cl),  s  G  R.  We  say  that 
u  =  (uo,Ui)  G  H1~s{ Cl)  satisfies  (13)  in  weak  sense ,  or  that  it  is  a  weak  solution  of 
the  Dirichlet  problem  (13),  if 

(16)  B(u,  v)  =  -(/,  v)n  +  ( g ,  duv)dn, 
for  all  v  G  H1+s(Cl). 

Remark  3.2.  If  u  =  (uq,ui)  G  H1~s( Cl)  is  a  weak  solution  of  the  Dirichlet  problem 
(13)  in  the  sense  of  the  above  definition  for  s  <  0,  then  it  is  easy  to  see  that  the 
first  component  uo  is  a  classical  solution  for  (13)  and  u-\  =  <9„Wo  on  dCl. 

Remark  3.3.  If  u  =  (uo,ui)  G  H1~s( Cl)  is  a  solution  of  (13)  in  weak  sense,  then  it 
is  unique  with  this  property. 

The  main  ingredient  we  will  use  for  proving  Theorem  0.3  is  the  following  “inf- 
sup”  lemma  (or  the  Babuska-Brezzi  condition)  [2,  3].  (This  result  was  used  for 
similar  purposes  in  [9]  in  order  to  deal  with  low  regularity  Neumann  data.) 

Theorem  3.4.  Let  X  and  Y  be  reflexive  Banach  spaces  with  norms  ||  •  ||x  and 
j|  •  || v-  Also,  let  B\  :  X  x  Y  — >  C  be  a  bilinear  form.  Assume  that 

(a)  B\  is  continuous; 

(b)  There  exists  7  >  0  such  that  inf||.E|ij^1  supm^n^  \B\{x,y)\  >  7; 

(c)  sup||a,||x<1  \B1(x,y)\  >  0  whenever  0. 

Then  for  any  continuous  functional  F  :Y  — >  C  there  exists  a  unique  x  G  X  such 
that  F(y)  =  Bi(x,y),  for  all  y  G  Y.  Moreover,  we  have  ||x||  <  ||F||/7. 

We  are  ready  now  to  prove  the  main  result  of  this  section,  that  is,  Theorem  0.3. 

Proof.  It  is  easy  to  see  that  B  is  continuous  from  its  definition  and  the  definition  of 
negative  order  Sobolev  spaces.  Therefore,  Condition  (a)  in  Theorem  3.4  is  satisfied. 
Let  u  =  (u0,ui)  G  £T1-S(fl)  be  such  that  ||u||  :=  ||uo||Ht—(n)  +  ||«i||i?-t/a— (an)  =  !• 
Since  u  also  belongs  to  (H~1+s(Cl)  ®  H1^2+s (dCl))* ,  which  is  the  dual  of  7L1-S(fl), 
there  exists  (v,v\)  G  H~1+s(Cl)  ®  Hl/2+s(dCl)  with  ||(u,ui)||  :=  ||w||a-i+»(n)  + 
||ui||ffi/2+S(0n)  =  1  such  that 

(17)  (u,  (v,  vi))  :=  (u0,  v)n  +  (ui,  tfi)an  >  1/2. 

Let  V  G  H1+s(Cl)  D  H1/2+s(dCl)  be  the  unique  solution  of  the  inhomogeneous 
Dirichlet  problem 


(18) 


AV  =  v  in  Cl, 

V  =  V\  on  dCl. 


Then, 

(19)  B{u,  V)  =  (u0,  AF)fi  +  (ui,  V)an  =  (w0,  u)n  +  (mi,  ^i)an  >  1/2, 
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and  this  inequality  implies  that  Condition  (6)  in  Theorem  3.4  is  also  satisfied. 

Finaly,  let  us  check  Condition  (c)  in  Theorem  3.4.  Let  v  £  H1+s( ft)  such  that 
B(u,v )  =  0,  for  all  u  =  (110,111)  £  H1~s( fi).  Then,  v  must  satisfy  the  Dirichlet 
problem 


(20) 


Av  =  0  in  il. 

v  =  0  on  dCl, 


whose  unique  solution  is  v  =  0.  This  shows  that  Condition  (c)  in  Theorem  3.4  is 
satisfied  as  well.  The  conclusion  of  the  Theorem  0.3  follows  if  we  take  the  functional 
F  in  Theorem  3.4  to  be  F(v)  :=  -(/,  v)q  +  ( g ,  duv)dn-  □ 


Part  2.  GFEM  Approximation  Spaces 

4.  The  Generalized  Finite  Element  Method 

Our  goal  is  to  construct  a  sequence  \i  =  1,2,...,  of  Generalized  Finite 
Element  spaces  that  satisfy  the  two  assumptions  of  the  previous  section.  To  this 
end,  we  shall  introduce  a  sequence  of  Generalized  Finite  Element  spaces  that  satisfy 
certain  conditions  (Conditions  A(/i;j),  B,  C,  and  D).  In  the  following  sections  we 
shall  prove  that  these  conditions  imply  Assumptions  1  and  2. 

We  begin  by  recalling  a  few  basic  facts  about  the  Generalized  Finite  Element 
Method  [4,  8,  26].  This  method  is  quite  convenient  when  one  needs  test  or  trial 
spaces  with  high  regularity. 

4.1.  Basic  facts.  Let  k  £  Z+.  We  shall  denote  as  usual 

Mw*.°°(n)  :=  max  ||<9au||Loo(n),  ||u||^,oo(n)  :=  max  ||aQu||Loc(n), 

|  a  | = Ac  |  o:  |  <  A; 

:=  {u,||u||w,*.«,(n)  <  00},  and  || Vw||^,oo(n)  :=  T,j\\dM\w°’°°(n)-  In 
particular,  |«|wo,~(a)  =  IM|iv°-“(n)  =  IM|l°°(J2)- 

We  shall  need  the  following  slight  generalization  of  a  definition  from  [8,  26]: 

Definition  4.1.  Let  C  R"  be  an  open  set  and  {uJj}jLi  be  an  open  cover  of  0 
such  that  any  x  €  O  belongs  to  at  most  n  of  the  sets  u>j.  Also,  let  be  a  partition 
of  unity  consisting  of  W/m,oo(0)  functions  and  subordinated  to  the  covering  {u>j} 
(i.e.,  supp  (f>j  C  tUj).  If 

(21)  ||9a^|Uoc(a)  <  C'fc/(dianiWj)fc,  k  =  |a|  <  m, 

for  any  j  =  1, . . . ,  N,  then  { <fij }  is  called  a  ( k ,  Go,  Gi, . . . ,  Gm)  partition  of  unity. 

Assume  also  that  we  are  given  linear  subspaces  C  j  =  1,2 , ...  ,7V. 

The  spaces  Tj  will  be  called  local  approximation  spaces  and  will  be  used  to  define 
the  space 

N 

(22)  S  =  Sgfem  ■=  {  Y,  <I>M,  v3  e  tfj}  c  Hm(n), 

3- 1 

which  will  be  called  the  GFEM-space.  The  set  {oJj,  (j)j.  Tj}  will  be  called  the  set  of 
data  defining  the  GFEM-space  S.  A  basic  approximation  property  of  the  GFEM- 
spaces  is  the  following  Theorem  from  [8]. 
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Theorem  4.2  (Babuska-Melenk).  We  shall  use  the  notations  and  definitions  of 
Definition  f.l  and  after.  Let  {<fj}  be  a  (k,Cq,Ci)  partition  of  unity.  Also,  let 
Vj  £  4'j  C  uap  :=  'Yhj4>jvj  £  <5.  and  dj  =  diamwj,  the  diameter  of  uij. 

Then 


(23) 


||V(u- wap)|||2(n)  <2 


u-uap |H2(n)  <  kC,^^||w-^|||2(w.)  and 

3 

cl  h-vAhw 


K)2 


C'ol|V(«-t;j)||i2(a 


<)  ■ 


4.2.  Conditions  on  GFEM  data  defining  5M.  Recall  that  uj  is  star-shaped  with 
respect  to  uj*  C  uj  if,  for  every  x  £  uj  and  every  y  £  lo*,  the  segment  with  end 
points  x  and  y  is  completely  contained  in  to.  Also,  recall  that  we  have  fixed  an 
integer  m  that  plays  the  role  of  the  order  of  approximation.  Let  {uij,  'hjlyLi 
be  a  single,  fixed  data  defining  a  GFEM-space  S,  as  in  the  previous  subsection, 
and  let  E  :=  {u>j,  (j>j,  fkj,  uj*},  where  u>j  is  star-shaped  with  respect  to  uj*  C  ujj.  We 
shall  need,  in  fact,  to  consider  a  sequence  of  such  data 


(24) 


E^  —  {o 


3  Jj=l’ 


p£N, 


defining  GFEM-spaces  S M 


(25)  :=  {  E  <%vi  ’  VJ  G  }  C  i? m  (0) , 

i=i 

such  that  there  exist  constants  A,  C, ,  cr,  and  k  and  a  sequence  h M  — >  0,  as  /i  — »  oo, 
for  which  Ep  satisfies  Conditions  A (/iM),  B,  C,  and  D  below  for  /t  €  N.  The  sequence 
hlt  gives  the  “typical  size”  of  the  elements  defining  S/t ,  as  in  the  first  part. 

Condition  A (/iM).  We  ftove  that  LI  =  uj1;'  and  for  each  j  =  1,2, .. . ,  7V/4 ,  the 

set  LOj  is  open  of  diameter  dj  <  <  1  and  uj C  is  an  open  ball  of  diameter 

>  adj  such  that  uj is  star-shaped  with  respect  to  uj*T  . 

Notice  that  we  only  assume  the  open  covering  {ujj}  to  be  nondegenerate,  a 
weaker  condition  than  quasi-uniformity  (see  [16],  Section  4.4,  for  definitions  and 
more  information  on  these  notions). 

Condition  B.  The  family  is  a  {n,  Co,  C\, . . . ,  Cm)  partition  of  unity. 

The  following  condition  defines  the  local  approximation  spaces  llC .  To  formulate 

this  condition,  let  us  choose  Xj  £  ujj  n  (XI,  if  the  intersection  is  not  empty.  We 
can  assume  that  linear  coordinates  have  been  chosen  such  that  Xj  =  0  and  the 
tangent  space  to  dLl  at  Xj  is  {xn  =  0}  =  R”_1.  For  hlt  small,  we  can  assume 
that  ujj  fl  dLl  is  contained  in  the  graph  of  a  smooth  function  g ^  :  Rn_1  — »  R.  If 
x  =  (x\,X2,  ■  ■  ■  ,xn)  £  R",  then  we  shall  denote  x'  =  {x\,X2,  ■  ■  ■  rxn-i)  G  Rn_1,  so 
that  x  =  (x',  xn).  Let  q1*  :  R™”1  — *  R  be  a  polynomial  of  order  m  such  that 


(26)  \9^(x,)-q^x,)\<C(dl;r+1  and 

| Vgf(x')  -  S7qj(x')\  <  C(d^)m  for  all  (x',x„)  £  uj 
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This  condition  is  satisfied,  for  instance,  if  dagh( 0)  =  daq^{ 0),  for  all  |a|  <  to.  In 
this  case,  the  m-degree  polynomial  q h  :  Kn_1  — >  K  is  uniquely  defined  by  the  afore 
mentioned  requirement. 

Next,  denote  by  q^  :  R"  — » IR”  the  bijective  map 

(27)  q?(x)  =  q?(x',xn)  =  (xr  ,xn  +  q^(x')). 

Let  us  denote  by  Vk  the  space  of  polynomials  of  order  at  most  k  in  n  variables. 

Condition  C.  We  have  A/h  =  Vm  if  D  dLl  =  0  and,  otherwise, 

=  {p  °  (qf  )“\  P  e  Vm,  such  that  p(x',  0)  =  0}, 

where  q h  are  polynomials  satisfying  Equation  (26)  with  a  constant  C  independent 
of  j  and  p. 

An  equivalent  form  of  the  condition  “p  £  Vm,  p{x' ,  0)  =  0”  is  “p  =  xnp\, 
Pi  G  Vm- 1,”  because  any  polynomial  vanishing  on  the  hyperplane  {xn  =  0}  is 
a  multiple  of  xn.  Since  {ql^)~1{x'  ,xn)  =  (x',xn  —  qh(x')),  we  obtain  p(x',  xn)  = 
(xn  ~  q$(x'))pi  o  (eft)"1. 

Condition  D.  We  have  (j>j  =  1  on  for  all  j  =  1, . . . ,  iVM  for  which  ^  0. 

The  constants  Cj,  a,  and  k  will  be  called  structural  constants.  Note  that  we 
must  have  N M  — >  oo  as  p  — >  oo. 

The  above  assumptions  are  slightly  weaker  than  the  ones  introduced  in  [9].  For 
instance,  Condition  C  implies  the  following  propety  (which  is  similar  to  Condition 
C  in  [9]) 

For  any  w  £  A>j,  any  0  <  l  <  m  +  1,  and  any  ball  uj*  C  toh  of  diameter  >  adk . 

(28)  IMI //'(<)  -  C\\wWh1(u>-). 

For  further  applications,  we  shall  also  need  a  variant  of  the  spaces  S,,  in  which  no 
boundary  conditions  are  imposed.  Recall  the  functions  q h  used  to  define  the  spaces 
Tj  .  Let  4 I1.  =  Afk  if  Uj  does  not  touch  the  boundary  dLl  and  Afj  =  {po  (qj)^1,  p  G 
Vm}  otherwise  (the  difference  is  that  we  no  longer  require  p  to  vanish  when  xn  =  0). 
We  then  define 

(29)  ^  ;=  {  E  <%vi  ’  vi  e  }  C  Hm (Q) . 

i= i 

We  shall  also  need  the  following  standard  lemma,  a  proof  of  which,  for  s  £  Z+, 
can  be  found  in  [9].  For  s  >  0  it  is  proved  by  interpolation. 

Lemma  4.3.  Let  ipj  be  measurable  functions  defined  on  an  open  set  W  and  s  >  0. 
Assume  that  there  exists  an  integer  k  such  that  a  point  x  £  W  can  belong  to  no 
more  than  k  of  the  sets  supp (ipj).  Let  f  =  Then  there  exists  a  constant 

C  >  0,  depending  only  on  k,  such  that  ||/|||/a(^)  <  C  W^iWh^w)- 

Recall  that  dj  denotes  the  diameter  of  (x'J- .  Let  us  observe  that  Condition  A(/im) 
implies  the  following  inverse  inequality. 
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Lemma  4.4.  There  exists  C  >  0,  depending  only  on  a,  such  that 
(30)  lbll*.K)  <  C(d^s\\p\\HrK), 

for  all  0  <  r  <  s  <  m,  all  j,  all  p,  and  all  polynomials  p  of  order  m. 

Proof.  The  proof  of  this  lemma  is  inspired  from  the  proof  of  (4.5.3)  Lemma  of  [16]. 
Consider  p  and  1  <  j  <  arbitrary,  but  fixed  for  the  moment.  Let 


where  x k  is  the  center  of  the  ball  to*h . 

If  p  G  Vm  is  a  polynomial  of  order  m,  then  p  is  defined  by  p(x)  :=  p(djX  +  x^ ) 
for  all  x.  Observe  that  the  set  Vm  :=  {p  :  p  £  Vm}  is  nothing  but  the  set  of  all 
m-degree  polynomials  in  x.  Clearly, 

(31)  \p\h*(q»)  =  {dj)k~n/2\p\HHu>f)’  for  0  <  k  <  m. 

We  first  prove  (30)  for  the  case  r  =  0.  Since  Vm  is  finite  dimensional,  we  have 
by  the  equivalence  of  norms  on  the  unit  ball  B( 0, 1)  that 

(32)  l|p||iffc(B(o,i))  <  C'IIpIU2(b(o,i)),  for  any  0  <  k  <  m, 

where  C  >  0  is  a  constant  that  does  not  depend  on  k,  j,  and  p.  From  Condition 
A (h/f),  we  obtain  that 

(33)  ||p||b2(b(o,i))  <  CWpWl^^) 

where  C  >  0  depends  only  on  the  structural  constant  a.  From  (32)  and  (33),  it  is 
clear  that 

IIpII <  C\\p\\L^)  \/p  G  Vm, 

where  C  >  0  depends  only  on  a.  Therefore,  (31)  implies 

|p|BH^)(di)fc_n/2  '5  C\\p\\L2(u,f)(dj)~n/2  for  0  <  k  <  s, 
from  which  we  deduce  that 

\p\h^)  <  C(dj)~k\\p\\L2(u»)  for  0  <  k  <  s. 

Since  dk  <  <  1,  we  have 

(34)  \\p\\h°(w»)  <  C(dj)~s\\p\\L2(u»), 
which  is  just  (30)  for  r  =  0. 

Let  us  now  analyse  the  general  case  0  <  r  <  s  <  m.  For  |a|  =  k,  with  s  —  r  < 
k  <  s,  Dap  =  DdD^p  for  \/3\  =  s  —  r  and  |q|  =  k  +  i —  s.  Therefore, 

II-d“pIIl2(uA‘)  <  WpWh.-^ 

<C[$y-'\\D'<p\\„{u>)  (by  (34)) 

<  c{qy-'\P\H*+r-.{u»). 

Since 

\p\Hk(ujf)  ;=  X!  \\Dap\\L^), 

\ot\—k 

we  obtain  that 

\p\h^)  <  C(d^)r~s\p\Hk+r-s(^)  for  s  —  r  <  k  <  s. 
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This  implies  that 

(35)  \p\h*(u,»)  <  C(djY~s\\p\\Hr(w»)  for  s-r<k<s. 

From  (34),  we  also  have 

(36)  lbll^K)  <  c{qr-'\\p\\L,{u!;)  < 

Combining  (35)  and  (36)  gives  (30)  and  this  ends  the  proof  of  the  lemma.  □ 


5.  Properties  of  the  spaces 

In  this  section,  we  establish  some  properties  of  the  GFEM  spaces  S^,  p  G  N, 
defined  in  Equation  (25)  using  the  data  =  {u)j,  ^ ,  'f'j,  w*M}j=i  satisfying  con¬ 
ditions  A (hY),  B,  C,  and  D  introduced  in  the  previous  section  for  —>  0.  The 
main  result  is  that  the  sequence  .S'M  satisfies  Assumptions  1  and  2  of  the  first  sec¬ 
tion.  Also,  we  prove  that  there  is  a  sequence  Gk  G  Sk  of  approximate  extensions 
of  g  which  satisfies  Assumption  3  in  the  case  of  the  non-homogeneous  Dirichlet 
boundary- value  problem  (13). 

Hereafter,  for  simplicity,  we  will  omit  the  index  p  whenever  its  appearance  is 
implicit. 

Let  us  fix  j  such  that  ZJ jP\dQ  is  not  empty.  Recall  the  functions  g-j ,  q3  :  Rn_1  — >  R 
defined  in  the  previous  section.  So,  for  h  small,  u}jC\d£l  is  contained  in  {(x1 ,  gj(x'))}, 
the  graph  of  the  smooth  function  gj  :  R"_1  — >  R  (this  may  require  a  preliminary 
rotation,  which  is  not  included  in  the  notation,  however,  for  the  sake  of  simplicity). 
Let  qj  :  R"  — »  R"  be  the  bijective  map  defined  by  Equation  (27).  Similarly,  let 

(37)  gj(x)  =gj(x',xn )  =  (a/,  xn  +  gj(x')). 

Then  gj  maps  R"_1  to  a  surface  containing  ujj  fl  dfl.  We  have  gj1{x)  =  (x',xn  — 
9j{x'))  and  qj1(x)  =  (x',x„  -  qj{x')). 

We  shall  need  the  following  estimate. 

Lemma  5.1.  For  any  polynomial  p  of  order  m,  we  have 

lb°  gj1  -poqj1  IU2K)  <  Cd™+1\\p\\m(wj)  and 
\\p°~gjl  -poq^Whh^)  <  Cd™\\p\\Hi{bJ]), 

where  C  is  a  constant  independent  of  p,  p,  and  j. 

Proof.  By  Taylor’s  expansion  theorem  in  the  xn  variable,  we  have 


P°9j  \x',xn)  =  p(x' ,  xn  -  gj(x'))  =  p(x',xn )  -  gj(x,)dnp(x' ,  xn)  +  ■■■ 


+  (-1)' 


k\ 


J\k 


dnP(x',  Xn)  +  .  .  .  +  (-1)’ 


,9 j(x'y 


d™p(x',xn) 


and 


P  °  Qj  1{rf,xn)  =  p(x',xn  -  qj(x'))  =  p(x',xn)  -  qj(x')dnp{x'  ,xn)  +  ... 


+  MV 


k\ 


J\k 


-d%p(x',xn)  +  . . .  +  (-1)’ 


iQj(x')r 


-dnP(x',Xn). 
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Then, 


I P°9j  (x',xn)  -p°qj  {x',xn)\  =  I p(x',xn  -  gj(x'))  -p(x',xn  -  qj(x')) \ 

<  \9j{x')  -  q-j{x')  |  •  \dnp{x',xn)\  +  . . .  +  •  \dkp(x',  xn)\  +  ... 


gj(xT  -  qj(x'Y 
to! 


I  dnP(x',Xn)\. 


From  this  and  the  Cauchy-Schwartz  inequality,  we  obtain 

I P°9j1{x',xn)  -poq~1(x',xn)\2  =  | p(x',xn  -  gj(x'))  -  p{x',xn  -  qj{x'))\2 
<  m[{gj{x')  -  qj(x'))2dnp{ x' ,  x„ )2  +  . . .  +  ^  ^  ^  dnP(x xn)2  +  . . 


(fc!)2 

{gj{x'r-qAx')m? 

(to!)2 


d™p(x',xn)2}. 


Notice  that  \gj(x')\  <  dj,  for  all  ( x\xn )  €  u>j,  and  because  qj{x')  =  gj( x')  + 
0(d™+1),  for  all  ( x',xn )  €  Wj,  we  have 

(g^x'f-q^x'f)2  =  ig*(x')-(gj(x')+0(d?+1))k}2  <  Cdfm+k\  for  k  =  1,. . 
which  in  turn  implies  that 

I P°gJ1(x',xn)-poq~1(x',xn)\2  <  Cdfm+1)[dnp(x',xn)2  +  d2d2p(x',xn)2  +  ... 

+  d2^1'* dkp(x' ,xn)2  +  . . .  +  d™p(x' ,xn)2]. 


By  using  the  inverse  inequality  dk  | |f>|| <  C\\p\\Hi^.),  we  get 


\\p°9j  1  ~P°q3  1\\l2(wj)  =  /  I P°9j  1(x',xn)  -P°qj  1{x' ,Xn)\2dx'dxn 

Jujj 

<  Cdfm+1)  f  [dnp(x',  xn)2  +  d2d2p(x',  xn)2  +  . . .  +  d2(k~1}  dkp(x' ,  xn)2  +  ... 

Jujj 

+  d™p(x' ,  xn)2]dx'  dxn 

<  Cdfm+1)[\\p\\~Hi^^  +  •  ■  ■  +  1'>\\p\\2Hk(Uj )  +  ■  •  •  +  dfm 

<Cdfm+1)\\pfHHuj), 

and  this  completes  the  proof  of  \\pogJ1  -  p  o  qJx\\u2(Uj)  <  CdJ+1\\p\\Hi^.y 
The  proof  of  \\pogJ1  —  P°qJ1\\m(uj:j)  <  C'^lbll^q^)  is  reduced  to  the  previous 
inequality  as  follows.  First,  from  the  inverse  inequality  dj|b||#i(w.)  <  C*|| j?|| ^2 (c^ - ) ? 
we  obtain 

(38)  WpogJ1  -poq^h^cv,)  <  CdJ\\p\\Li(uJ:i). 

It  is  then  enough  to  show  that 

(39)  WdkipogJ1)  -  dk(poqT1)\\L2{u.)  <  Cd™\\p\\Hi(U]), 
for  all  k  =  1,2, ...  ,n. 
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The  case  k  =  n  is  easier,  so  we  shall  treat  only  the  case  when  1  <  /c  <  n  —  1.  A 
Taylor  expansion  with  respect  to  the  ccn-variable  gives 

dk(pog~1)(x',xn)  =  dk(p(x',xn  -  gj(x'))) 

=  (■ dkp)(x',xn  -  gj{x'))  -  dkgj{x')(dnp){ x',xn  -  gj(x')) 

and 

dk{p  °  ,xn)  =  dk(p(x',xn  - 

=  ( dkp)(x',xn  -  qj( x'))  -  dkqJ{x')(dnp){x' ,xn  -  qjfx')) 

Equation  (39)  then  follows  from  Equation  (38)  and  from  the  estimates  qj{x')  = 
gj(x')  +  0(dj+1 ),  dkqj( x')  =  dkgj(x’)  +  0(d J1)  and  \gj{x')\  <  dj  for  (x',xn)  G  uj 
(see  Equation  (26)  and  Condition  C).  □ 

Remark  5.2.  Let  us  observe  that  Condition  A (/iM)  was  used  implicitly  in  the  proof 
of  Lemma  5.1  when  we  used  the  inverse  estimates  d^1\\p\\Hk^UJj')  <  C'||p||/ji(aiJ)- 

Remark  5.3.  If  (26)  is  replaced  by  the  more  restrictive  condition  \da(gj  —  qj)\  < 
Cd™+1~ l“l,  for  all  |a|  <  m+ 1,  then  the  result  of  the  above  lemma  can  be  extended 
as  follows:  For  any  polynomial  p  of  order  m,  we  have 

WpogJ1  -P°qj1\\H-{wj)  <  CdJ+l~s\\p\\Hi^.),  s  =  0, . . .  ,m+  1, 
where  C  is  a  constant  independent  of  p,  /./,  j,  and  s. 

Corollary  5.4.  Let  p  G  Vm,  then 

II <Po{p°9jl  -p°971)IIhik)  ^  Cd?\\p\\Hi(Uj)- 

If  p  G  Vm  also  vanishes  on  {xn  =  0}  then  we  have 
\\<Pj(P°  ?71)IU2(9«)  - 

Here  C  is  a  constant  independent  of  p,  p,  and  j. 

Proof.  Using  Lemma  5.1  and  Assumption  B,  we  obtain 

II faipogj1  -poqf1) Whhwj)  <  \pogJ1  - p° qj'Wm^) 

+  (xuildihh-^llpo-gj1  -poq^hn^) 

<  Cd?\\p\\HHuj)  +  Cdj^1  \\p\\HHuj)  <  Cd?\\p\\HHwj). 

The  last  part  follows  from  the  first  part  of  this  corollary,  which  we  have  already 
proved,  and  from  the  fact  that  4>j(p°  gj1)  =  0  on  dLl.  Indeed, 

ll<A?(p°  qJ1)\\L2(dn)  =  Hj{p  °  gj1  -  p  °  qj^h^dn) 

<  C\\ (fjipogj1  -peqJ^Wmin)  =  CWfaijpogJ1  -pog71)||Hi(w.) 

<Cdf\\p\\m{uj)<Cd^\\p\\m{u)V 

The  proof  is  now  complete.  □ 

We  are  ready  now  to  prove  that  Assumption  1  is  satisfied  by  the  sequence  of 
GFEM- spaces  Stl  introduced  in  Subsection  4.2. 

Proposition  5.5.  Let  S ^  be  the  sequence  of  GFEM-spaces  defined  by  data 
(Equation  (24),)  satisfying  conditions  Alh^),  B,  C,  and  D.  Then  the  sequence  S M 
satisfies  Assumption  1. 
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Proof.  Let  Wj  £i&j  and  w  =  Y2  (t>jwj  £  Since  we  are  interested  in  evaluating  w 
at  dtt,  we  can  assume  that  only  the  terms  corresponding  to  j  for  which  uJj  fl  Oil  0 
appear  in  the  sum.  Then  wj  =  Pj  °  qj  ,  for  some  polynomials  pj  £  Vrn  vanishing 
on  {xn  =0}.  Hence  Lemma  4.3  and  Corollary  5.4  give 


Mlifi/^en)  <  ll^u,jllia(en)  =  °Qj  1)lli=(an) 

3  3 

<  Ell^ll^K)  ^  C^EIKHW)- 


=  w 


By  Condition  D,  \\wj\\m(u*) 

\H\i  Hon)  <  chim\\wrH1(u w;)  <  ch^r\\wrHi{n). 

Assumption  1  is  hence  satisfied  by  taking  square  roots. 


H1(Uwj 
2 


y  Therefore, 

'  2mii  II 2 


□ 


Remark  5.6.  Condition  D  is  only  needed  in  the  proof  of  Proposition  5.5.  Although 
one  can  prove  that 

(40)  E  Ill'll —  KKIIfl-qn) 

3 

(by  using  norm  equivalence  in  finite  dimensional  spaces),  one  can  not  bypass  Con¬ 
dition  D  because  the  constant  C  in  (40)  depends  on  p.  To  remove  this  dependence, 
one  would  have  to  impose  additional  and/or  different  conditions  on  the  partition 
of  unity. 


The  proof  that  the  sequence  5/,  also  satisfies  Assumption  2  is  also  based  on 
the  above  lemma  and  on  the  following  result.  Recall  that  the  local  approximation 
spaces  T ,  and  T/  were  defined  in  Subsection  4.2. 

J  3 


Lemma  5.7.  Let  u  £  Hm+1(uij).  Then  there  exists  a  polynomial  w  £  such 
that  \\u  -  w\\Hi{uj)  <  Cd™\\u\\Hm+i{Uj)  and  \[u  -  w\\L2{Uj)  <  Cd™+1\\u\\Hm+ i{ujj) 
for  a  constant  C  independent  of  u,  p,  and  j.  If  u  =  0  on  uij  fl  dft,  then  we  can 
choose  w  £  'L/ . 


Proof.  We  are  especially  interested  in  the  case  when  u  =  0  on  ZJj  fl  dCl,  so  we  shall 
deal  with  this  case  in  detail.  The  other  one  is  proved  in  exactly  the  same  way. 

Let  us  consider  v  =  u  o  pj.  Since  gj  maps  R””1  =  {xn  =  0}  to  a  surface 
containing  UJj  fl  dfl,  we  obtain  that  v  =  0  on  R"_1.  For  small  enough,  we  can 
assume  that  g^Kj)  lies  on  one  side  of  R™_1.  Let  U  be  the  union  of  the  closure  of 
5j1(wj)  and  of  its  symmetric  subset  with  respect  to  R”_1.  Define  v\  £  H 1(U)  to 
be  the  odd  extension  of  v  (odd  with  respect  to  the  reflection  about  the  subspace 
R"_1).  Let  pi  be  the  projection  of  iq  onto  the  subspace  Vm  of  polynomials  of 
degree  m  on  U .  This  projection  maps  even  functions  to  even  functions  and  odd 
functions  to  odd  functions.  Hence  p\  is  also  odd.  In  particular,  p\  =  0  on  Rn_1. 
We  also  know  from  standard  approximation  results  [16]  that 

IK  _  Pill  J?1  (17)  <  Cdjl\\v1\\Hm+i(Uy 

Then 


\\u-Pi°9j  < 


^  II  ^ i  y i.  i 
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Let  w  =  pi  o  q-  1 .  The  lemma  follows  from 

ll«- <  \\u-pi  °5“1||h1(^)  +  ||pi  ogj1  -Pi 

<  Cdp\\u\\Hm+i{Uj)  +  Cd?\\Pl\\Hiiuj)  <  Cd?\\u\\Hm+Huj), 

where  we  have  also  used  Lemma  5.1. 

To  prove  the  relation  ||u  —  w||l2(w3)  <  Cd™+1||w||#m+i(a,i),  we  first  notice  that 
Poincare’s  inequality  gives 

IK  -  Pi\\l2(u)  <  C'd,- ||«i  Pi Wh1  (u)  <  C'd™+1||iii||ffm+i(C/). 

The  rest  is  exactly  the  same.  □ 

We  are  ready  now  to  prove  Assumption  2.  See  [4],  section  6.1,  and  [9]  for  related 
results. 

Proposition  5.8.  The  sequence  of  GFEM  spaces  S M  satisfies  Assumption  2. 

Proof.  We  proceed  as  in  [9],  Theorem  3.2.  Let  u  €  Hm+1(PL).  If  UJ.j  does  not 
intersect  <9f l,  we  define  Wj  £  \Pj  =  Vm  to  be  the  orthogonal  projection  of  u  onto  Vm 
in  IL1(wj).  Otherwise,  we  define  Wj  £  Afj  using  Lemma  5.7.  Then  let  w  =  4>jwj- 
By  using  Lemma  5.7,  the  definition  of  the  local  approximation  spaces  ’Lj  (Condition 
C),  and  the  bounds  on  ||V</>j||L°°  (Condition  B),  we  obtain 

I  u  -  w  |  Hi  (n)  =  |  <t>j  (u  -  wi )  I  m  (Q) 

3 

3 

<  E  (CdT KII^+h^)  +  cdj'd^M^1^)  <  CKh™\\u\\Hm+1(u). 
j 

This  completes  the  result.  □ 

Next,  we  will  be  looking  for  a  sequence  Gk  of  approximate  extensions  of  g  in  Sk- 
Recall  that  the  spaces  Sk  D  Sk  were  defined  in  Equation  (29)  and  are  variants  of 
the  spaces  Sk  that  are  not  required  to  satisfy,  even  approximately,  the  boundary 
conditions.  The  construction  of  such  a  sequence  Gk  of  approximate  extension  as 
well  as  the  analysis  of  the  resulting  method  are  the  main  results  of  this  subsection. 
Other  methods  for  constructing  Gk  are  certainly  possible. 

We  now  check  that  it  is  possible  to  choose  Gk  €  Sk  satisfying  Assumption  3.  We 
follow  the  method  in  [3]. 

Proposition  5.9.  There  exist  continuous  linear  maps  Ik  :  Hm+1( ft)  — >  Sk,  such 
that 

(41)  \u  —  Ik(u)\Hr(Q)  <  Ch™  +  1  r|H|_F/m.+i(Q), 

for  r  =  0  and  r  =  1. 

Proof.  For  u  £  Hm+1(fl)  and  j  fixed,  let  v  =  u  o  gj.  The  Taylor  polynomial  of 
degree  m  of  v  averaged  over  <)-1(wj)  is  given  by 

Pj{x)  :=  Qfv(x)  =  j 

Jg  1(uj) 


(42) 


Qy,v,n{^^j  (?/)  dy, 
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where 

'  v^(y) 

Qy,v,n{x )  =  v(y)+y^jdiv(y)(xi-yi)  +  . . .+  ^  - — j -(x-y)a,  a\  =  a  i!...aj, 
i= 1  |a|=m 

is  the  Taylor  polynomial  of  v  at  y  of  degree  m  and  G  C^°{g~1{iOj))  is  a  function 
with  integral  1.  Then,  by  the  Bramble-Hilbert  Lemma,  we  have 

(43)  \v  -  Pj\Hs{s-i{Uj))  <  Ch™+1~s \v\Hm+i(§-i („.)),  for  all  0  <  s  <  m+  1. 

Consider  Wj  :=  Pj  o  gT1  G  Tj.  Let  w  :=  Yhj  (l>jwj-  Then, 


(44)  \u-w\2Hr{n)  <  C^2\4>j{u-w)\2Hr{n)  <  CY^\Mu~w)\H'(Ui) 


^  lrvi’oo(wJ) \u  wj\ 2Hr~i(ujj) 

j  *=o 

r 

—  C  \(l)j\'wi’°°(u>j)i\U  ~  Pj  °  9j  \2HT-i(uj)Jr\Pj°9j  —  Pj  °  Qj  I 

j  *=o 

By  changing  variables  and  (43),  we  obtain 

(45) 

\u  -  Pj  o  gj1  \2Hr-i{Uj)  =  \vo  gj1  -  Pj  o  ~gj1\1Hr-i{uj])  <  C\v  -  Pj 


<  ch2(m+t~r+i^  |»''2 


=  Ch 


2(m+l— r+i) 


<  2(m+l  — r+z)  |  1 2 

-  °'lfc  lUl Hm+1(wj)- 


Also,  from  Lemma  5.1  and  the  definition  (42)  of  Pj ,  we  have 

(46)  | Pj  o  j"‘  -  Pj  o  <  CA?”+,-r+‘,||PJ||2H1(„.) 

From  (44),  (45),  (46),  and  Condition  B,  it  follows  that 

(47) 


i«-»iW)ScEEftC[*r+1-r+,)i<*i?,-+.(-,)+'‘r"-r™w?,.,»J,] 


-  2(m+l-r+i)  m  m2 


j  i= 0 


<  ChHm+1~r)  J2  IMI <  Cnhfm+1-r\\u\\l 


for  all  0  <  r  <  m  +  1. 

Define  Ik(u)  :=  uj.  Clearly  /*.  is  a  linear  map  from  Hm+1(Q)  to  Sj-  which  satisfies 
(41).  This  ends  the  proof  of  the  proposition.  □ 


Remark  5.10.  If  we  assume  the  stronger  condition  stated  in  Remark  5.3  on  the 
m-degree  polynomial  qj,  then  the  conclusion  of  Proposition  5.9  is  valid  for  0  <  r  < 
in  +  1  (the  proof  being  exactly  the  same) . 

From  this  proposition  we  obtain  right  away  the  Assumption  3. 

Proposition  5.11.  For  any  g  G  Hm+1/2(dQ)  there  exists  a  sequence  Gk  G  Sk 
satisfying  Assumption  3. 
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Proof.  Indeed,  let  us  chose  G  £  Hm+1(fi )  that  extends  g  to  the  interior  and 
satisfies  ||G||#m+i(n)  <  C\\g\\Hm+i/2^dQip  with  C  independent  of  g.  Then  choose 
Gk  =  Ik{G ),  with  Ik  as  in  Proposition  5.9.  □ 

6.  Interior  numerical  approximation 

In  this  section,  we  construct  a  sequence  of  approximations  Uk  G  Sk  of  the  solution 
u  of  the  distributional  boundary  value  problem  Aw  =  0  in  Q,  u  =  g  on  dfl,  with 
H1/2  s(<9f 1),  and  we  prove  interior  estimates  for  the  error  u  —  Uk ■  The  sequence 
of  spaces  Sk  is  the  sequence  of  GFEM  spaces  constructed  in  Section  4  and  hence  it 
satisfies  the  Assumptions  2  and  3,  by  the  results  of  Section  5.  In  particular,  Sk  C  Sk- 
We  need  to  consider  GFEM  spaces  instead  of  the  more  general  framework  of  the 
first  part  because  we  need  the  interior  approximation  result  of  [9]  recalled  below. 

In  this  section,  s  £  Z+  is  fixed.  Our  results  mirror  the  ones  in  [9],  where  the 
Neumann  problem  was  considered.  The  approach  is  different  however,  in  part  be¬ 
cause  it  is  more  involved  to  construct  finite  element  approximations  of  the  solution 
u  in  the  case  of  the  Dirichlet  boundary  conditions. 

Our  approach  is  to  first  approximate  g  with  a  sequence  Gk  of  functions.  Then 
each  of  the  approximate  equations  Aw*,  =  0,  w*,  =  Gk  at  the  boundary,  is  solved 
approximately  using  the  results  of  the  previous  sections.  This  yields,  for  any  k,  a 
sequence  (w*, j;i  with  (w*.)/t  £  S The  sequence  of  approximations  is  then  u M  := 
Our  approach  is  thus  similar  to  that  of  Section  2. 

Recall  the  spaces  S M  and  S^  introduced  in  Assumptions  1,  2,  and  3.  Let  us 
also  recall  that  S^  C  5M.  The  definition  of  the  space  S^  is  slightly  different  from 
the  one  in  [9];  however,  the  difference  is  only  in  the  local  approximation  spaces  at 
the  boundary,  and  hence  this  does  not  affect  the  spaces  S<( fi)  :=  S M  nCc°°(!l). 
Therefore,  it  can  be  shown  that  Theorem  3.12  from  [9]  is  still  valid.  Namely,  we 
have  the  following: 

Theorem  6.1.  Let  A  <g  B  C  be  open  subsets.  Then  there  exists  C  >  0  with  the 
following  property.  If  u  £  H1(  ft)  and  u M  £  S M  are  such  that  B(u  —  wM,x)  =  0  for 
all  x  G  S<( fl),  then  for  h „  small  enough, 

Wu-uJh^a)  <C[  inf  \\u-x\\m(B)  +  \\u-ull\\H-m{B)y 
The  constant  C  depends  only  on  A,  B,  and  the  structural  constants,  but  not  on 

/i  €  . 

This  result  is  the  version  for  the  Generalized  Finite  Element  of  a  basic  result 
by  Nitsche  and  Schatz  [29,  Theorem  5.1].  See  also  [31,  Theorem  9.2].  The  above 
result  is  the  reason  why  we  work  in  this  section  in  the  framework  of  the  Generalized 
Finite  Element  Method  and  not  in  the  abstract  setting  of  the  first  part. 

Recall  that  s  is  fixed  in  this  section.  We  otherwise  use  the  notation  of  the 
previous  sections.  We  shall  need  the  following  property  of  the  spaces  Sj. 

•  The  low  regularity  approximate  extension  property:  There  exists  a  constant 
C  >  0  such  that,  for  any  g  £  H1l'1~s(fd Q),  we  can  find  a  sequence  Gj  £  Sj 
satisfying  HG^n  -  <  ch]\\g\\Hi/z-s(dn)  for  all  0  <  7  <  m 

and  ||Gj||/fi(n)  <  ChJ!’\\g\\Hi/2-^any 
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Lemma  6.2.  Assumptions  2  and  3  imply  the  low  regidarity  approximate  exten¬ 
sion  property.  In  particular,  the  spaces  Sj  satisfy  the  low  regidarity  approximate 
extension  property. 

Proof.  Let  g  G  H1/2~s(dLl).  As  in  [30],  we  shall  consider  for  each  H  >  0  functions 
G#  G  C°°(Ll)  such  that 

II 9  ~  GhWjjI/2-s—!^)  <  CH'r\\g\\Hl/2-s(dty 

II Gh- || #1(0)  ^  GiJ  s||fl,ll#1/2-s(ao)., 

with  constants  independent  of  g.  On  R”,  we  can  chose  G#  =  \h  *g,  with  xh(x)  = 
H~nx{x / H)  for  a  suitable  chosen  x  G  C[?0(]R11).  In  general,  this  procedure  can  be 
localized  as  in  [15]. 

Let  us  then  choose  H  =  hj  and  Gj  G  Sj  to  have  the  following  projection  property 

B(Gh  —  Gj,x)  =  0,  x  £  Sj,  and  f  (G#  -  Gj)dS(x)  =  0. 

Jon 

Then  ||Gj||#i(n)  <  ||G#||#i(n)  <  G/i“s||5||#i/2-s(an).  To  prove  our  result,  it  is 
enough  to  show  that  the  restrictions  to  the  boundary  df l  satisfy 

(48)  || Gj  -  G#||ffi/2-i(ao)  <  G/ij||G#||#i(Q), 

for  any  l  >  0  (for  our  result  we  then  take  l  =  s  +  7). 

We  shall  proceed  as  in  [3].  Let  w  G  Hl~1^2(d Ll)  with  JdQwdS(x)  =  0.  We  let 
W  G  Hl+1(Ll)  be  the  unique  solution  with  mean  zero  of  the  Neumann  problem 
AW  =  0,  dvW  =  w.  Then  ||IL||#i+i(f2)  <  G||u>||fl-!-i/2(af2)  and  hence  Assumption  2 
(approximation  property)  gives  that  there  exists  y  G  Sj  such  that  \W  —  x\h1(q.'\  < 
Chj\\W\\Hi+i(n)  <  Chlj\\w\\Hi-i/2(dny  We  obtain 


I  Gj  -  G#||#i/2-i(am  —  sup 

0 


(Gj  -  Gh ,  w)oo 


=  sup 

w\\Hl-1/2(dQ.)  WjiO 

<  C\Gj  -  Gh |#i(n)  sup  ^  ~  A;|gl(n)  <  Ch) ||G#||#i(a). 

0  Ww\\Hl-1/2(dn) 

In  view  of  the  fact  that  Jda(Gj  —  Gn)dS(x)  =  0,  this  proves  Equation  (48)  and 


B(Gj-GH,W-X) 


hence  completes  the  proof  of  the  Lemma. 


□ 


We  are  now  ready  to  prove  a  result  on  the  interior  approximation  properties  for 
boundary  value  problems  with  low  regularity  (i.e.,  distributional)  boundary  data. 

Theorem  6.3.  Let  S M  C  Sfl  be  our  sequences  of  GFEM  spaces.  Also,  let  g  G 
Hl/2~s(dLl),  1  <  s  <  m  — 1,  Ad  LI,  and  Gj  be  as  in  the  low  regularity  approximate 
extension  property.  We  define  ( Vj )M  to  be  the  discrete  solution  of  the  equation 
Ait  =  0 ,  u  =  Gj  at  the  boundary,  as  defined  in  Equation  (54)  and  :=  (t;M)M. 
Then 

llM  -  *vll m(A)  <  Ch™~s~l\\g\\Hi/2-s(dn). 

Proof.  We  now  proceed  as  in  Section  2.  Let  us  denote  by  Wk  the  solution  of 
(49)  —A Wk  =  A Gfc  G  L2(Ll)  in  Ll  and  Wk  =  0  on  dLl, 

Then  we  let  Vk  :=  Wk  +  G^,  which  will  satisfy  Avk  =  0  in  LI,  Vk  =  Gk  on  dll. 
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Lemma  6.2  shows  that  the  low  regularity  approximate  extension  property  is 
satisfied.  Theorem  0.3  and  the  low  regularity  approximate  extension  property  then 
give 

(50)  ||u  —  Ufe||ffi-S-T(n)  <  C'H <7  —  Gfc||Hi/2-s-7(aa)  <  Ch'l\\g\\Hi/2-S(dQjy 

Let  B  be  an  open  set  such  that  A  <s  B  <s  and  t  be  a  parameter.  Then  for 
any  harmonic  function  tp  G  C°°(fl),  there  exists  a  constant  C  such  that  ||^||_yt(s)  < 
C\\(f>\\H1-i!-i(n)-  Taking  <p  :=  u  —  Vk,  t  =  1,  and  using  also  Equation  (50),  we  obtain 

(51)  ||u  -  Wfe|| Hi(B)  <  C\\u  -  Vfe||j?i-«-T(n)  <  Chl\\g\\Hi/2-S(dny 
By  taking  <j>  =  Vk,  we  obtain 

(52)  Kllffm+1(B)  <  C'|Kll.ff1-«(fi)  <  C'l|Gfc||ffi/2-s(an)  <  C\\g\\Hi/2-s(dny 
Also,  let  us  denote  by  ( Wk )^  G  S M  the  discrete  solution  of  the  problem  (49)  Namely, 

(53)  B((wk)n,x)  =  (AGfc,x),  X  G  K- 
This  is  nothing  but  Equation  (2)  for  /  =  0.  Let 

(54)  («fc)M  :=  {wk)n  +  Gfe. 

Then  Theorem  1.6  gives 

(55)  IK  -  KVIk-dO)  <  ChP+l+1\\AGk\\HP-Hn)  <  Ch?+l+1\\Gk\\HP+Hn), 

for  p  +  l  +  1  <  m. 

Let  us  now  take  p  =  0  and  l  =  s  +  7  —  1,  which  satisfy  s  +  j  =  p+  l  +  l<m. 
Then 

(56)  K  -  <  Chs^\\Gk\\HHn)  <  Ch^h^\\g\\H,n-s(dny 

In  particular,  for  k  =  g,  we  obtain 

(57)  IK  -  K)mI li?1-°-T(n)  <  Ch'p\\g\\Hi/2-s(dn)- 

Since  B{y^  -  K)M,  x)  =  B(w M  X)  =  0  for  all  x  €  Theorem  6.1  then 

gives 

(58)  |K  -  (^)Jffi(A)  <  C[  inf  IK~XlK(B)  +  IK~KKIKm (B)) 

<  G(K  +  fyK)IIS,lliL1/2-a(,9Q)  <  G/l^||5'||J?l/2-»(9n) 
where  the  last  two  terms  were  estimated  using  Equations  (52)  and  (57). 

Equations  (51)  and  (57)  give  for  s  +  7  +  1  =  m  that 

\\u  -  K)m|K(a)  ^  llu  -  v»\\m(A)  +  IK  -  K)mI \hha)  <  Chl\\g\\Hi/2-s{m) 

The  proof  is  complete.  □ 

7.  Comments  and  further  problems 

In  spite  of  all  the  differences  in  assumptions  and  definitions  between  [11,  28] 
and  our  paper,  the  main  issue  seems  to  be  providing  simple  examples  of  spaces 
satisfying  the  various  assumptions  used  in  these  papers.  For  instance,  it  would  be 
interesting  to  provide  other  examples  of  spaces  Slt  satisfying  Assumptions  1  and 
2.  It  would  also  be  interesting  to  see  if  a  modification  of  the  uniform  partition  of 
unity  can  give,  by  restriction,  spaces  .S'M  satisfying  these  Assumptions.  Finally,  it 
would  be  important  to  integrate  our  results  with  the  issues  arising  from  numerical 
integration  and  to  provide  explicit  numerical  examples  testing  our  results.  In  fact, 
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related  numerical  tests  together  with  some  theoretical  results  can  be  found,  for 
example,  in  [5,  20,  30]. 

7.1.  Comments  on  the  approximate  boundary  conditions.  For  most  of  our 
results,  it  is  enough  to  assume 

||*VlU2(9fi)  ^  Ch™ for  any  v M  €  S». 

However,  one  then  has  to  replace  ||w||#i(n)  with  ||u||,H-2(n)  in  Lemmas  1.4  and  1.5. 
Also,  the  proof  of  Theorem  1.6  for  p  =  0  requires  the  full  strength  of  Assumption 
1.  The  case  p  =  0  is  the  one  needed  for  the  results  of  Section  6,  and  that  is  the 
only  real  reason  for  which  we  need  the  stronger  form  of  Assumption  1. 

A  related  problem  is  to  construct  other  examples  of  spaces  satisfying  the  interior 
estimates  used  in  Section  6  as  well  as  the  general  Assumptions  1,  2,  and  3.  Then 
Theorem  6.3  would  be  valid  for  these  spaces  as  well. 
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